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The transverse self- modulation of ultra-relativistic, long lepton bunches in high-density plasmas 
is explored through full-scale particle-in-cell simulations. We demonstrate that long SLAC-type 
electron and positron bunches can become strongly radially self-modulated over centimeter distances, 
leading to wake excitation in the blowout regime with accelerating fields in excess of 20 GV/m. We 
show that particles energy variations exceeding 10 GeV can occur in meter-long plasmas. We 
find that the self-modulation of positively and negatively charged bunches differ when the blowout 
is reached. Seeding the self-modulation instability mitigates the effect of the competing hosing 
instability. This work reveals that a proof-of-principle experiment to test the physics of bunch 
self-modulation can be performed with available lepton bunches and with existing experimental 
apparatus and diagnostics. 
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I. INTRODUCTION 

Plasma-based acceleration offers the possibility to ac- 
celerate electrons and positrons in linear accelerators 
to high energies in distances much shorter than those 
achievable with conventional radio-frequency accelerat- 
ing devices. Laser pulse (laser wakefield acceleration or 
LWFA pQ) or particle bunch (plasma wakefield acceler- 
ation or PWFA [2]) drivers can excite large amplitude 
plasma waves or wakes. Peak electric fields in excess of 
E = 96 Vnolcm- 3 ] [V/m] can be generated, and accel- 
erating gradients in excess of 100 GV/m can be achieved 
using plasma densities no in the 10 18 cm -3 range, typical 
of current experiments. Plasmas may thus contribute to 
a novel generation of compact accelerators, with impact 
on many scientific and technological applications [3]. 

Proton bunches are attractive for plasma wakefield ac- 
celeration because the total amount of energy they carry 
(more than 100 kJ for a single LHC bunch) and the en- 
ergy per particle (7 TeV at LHC) is much larger than that 
of the lepton bunches of a future linear collider (1.6 kJ 
for a typical 2 x 10 10 , 500 GeV lepton bunch). One 
can therefore envisage accelerating a single ~ 10 GeV 
incoming lepton bunch to the TeV energy scale in a sin- 
gle PWFA stage driven by a relativistic proton bunch. 
The initial proposal for a proton-driven PWFA (or PDP- 
WFA) [4 considered the case of a short 100/mi) pro- 
ton bunch driver. However, such short proton bunches 
are not currently available. Kumar et al. [5 recently 
suggested that the radial self-modulation (S-M) of a long 
proton bunch at the wavelength of the relativistic plasma 
wave can lead to the resonant excitation of large ampli- 
tude plasma wakefields. This is analogous to the self- 
modulated LWFA [6 that was used in early experiments 
when fs-long laser pulses were not available. The self- 
modulation of a particle bunch is completely analogous 



to that of a laser pulse (photon beam) [7]. Studying ex- 
perimentally the physics of the S-M of particle bunches 
that could be considered as drivers for large energy gain 
PWFA experiments is therefore interesting and impor- 
tant. 

The S-M is the result of a transverse two-stream in- 
stability. For this instability to develop for high energy- 
beams the bunch transverse dimensions must be on the 
order of, or smaller than the cold plasma collisionless skin 
depth (a r « c/u p ) to avoid possible transverse bunch fil- 
amentation [8]. Since the S-M is a convective instabil- 
ity, longer bunches, which encompass a higher number of 
plasma skin depths (a^ ^> c/uo p ) and longer propagation 
distances (z ^> k^ 1 ) lead to larger wakefield amplitudes. 

Here c is the speed of light, uj p = \/n§e 2 / (ra e eo) the 
electron plasma frequency, m e the electron mass, e the 
elementary charge, eo the vacuum dielectric constant, z 
the propagation distance, hp = cj p /(cv // 27o)(^e/^p) 1 ^ 2 
the betatron wavenumber with m p the bunch particles' 
mass, jorripC 2 the energy of each beam particle, 70 the 
relativistic factor, and a r and the bunch rms radius 
and length, respectively. 

In the simulations presented here, we consider the 
case of high-energy electron and positron bunches that 
are currently available for experiments. We find that 
with these ultra-relativistic (70 oc 40,000) and long 
((j£ oc 500 fim) lepton bunches available at SLAC FACET 
[9] the saturation of the S-M instability is reached over 
only a few centimeters of plasma at a density of 2.3 x 
10 1T cm -3 . This is due to the lower relativistic mass (25 
GeV) of leptons when compared to that of protons in 
available bunches (450 GeV, See Table [i]). These lepton 
beams can therefore be focused tightly and large plasma 
densities can be used (a r « c/uo p oc n e ), which lead 
to large transverse focusing and longitudinal accelerat- 
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ing wakefield amplitudes. In addition the betatron wave- 
length (2ir/kp) is much shorter so shorter propagation 
distances are needed. 



In 3D numerical simulations we find that, as was previ- 
ously noted [10] , hosing-like instability of the bunch that 
is long when compared to the plasma wavelength may 
be a serious limitation to its stable propagation over dis- 
tances longer than the saturation length of the S-M in- 
stability. The hosing also occurs at a longer wavelength 
than that of the wake. Similar behavior was observed 
in simulations of long laser pulses [TTJ [12] . In an accel- 
erator scheme this may also limit the energy gain by a 
short, externally injected witness bunch. However, these 
simulations also indicate that the seeding of the S-M in- 
stability can significantly reduce the growth of the hos- 
ing instability over the plasma lengths considered here, 
as was noted in the laser pulse case [TTJ [12] . We con- 
sider the case of an electron bunch with an asymmet- 
ric current profile with rapidly rising edge to seed the 
S-M instability and defer to another publication the de- 
tailed study of the competition between the two instabil- 
ities. Further simulations with the shaped bunch in 2D 
cylindrically symmetric geometry that precludes hosing- 
like instabilities were performed to scan parameters. At 
saturation of the instability, the blowout regime [13] is 
reached and accelerating fields in excess of 20 GV/m are 
excited. However, in this nonlinear regime the behav- 
ior of electron and positron bunches differ, in contrast 
with the linear regime for which their behavior is iden- 
tical. Differences in wakefield excitation by short, pos- 
itively and negatively charged bunches were described 
previously [14]. Because most of the volume of plasma 
accelerating structure in the nonlinear regime (ion col- 
umn) is defocusing for positrons, a majority of positrons 
are quickly and strongly defocused and do not participate 
in the wakefields excitation. With the stable propagation 
of the electron and positron bunches over a meter-scale 
plasma, peak energy gain and loss by the bunch particles 
in the multi-GeV range are observed. 

This paper is organized as follows. In Sec. [II] we point 
out the differences between proton and lepton bunch S- 
M, and describe a possible S-M experiment using lepton 
bunches currently available at SL AC-FACET. In Sec. Ill 
we describe the base line parameters for the 3D and for 
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the 2D cylindrically symmetric simulations. In Sec. IV 



present 3D simulation results illustrating the competition 
between the hosing-like and the S-M instabilities with 
symmetric (with long rise time) and asymmetric (with 
rapid rise time) current profiles. In Sec. |V| we show 2D 
cylindrically symmetric simulation results of of the self- 
modulation of long SLAC electron and positron bunches. 
In Sec.[yi[ we describe the S-M wakefields excitation over 
meter long plasmas. In Sec. |VIII| we refer to possible 
diagnostics to observe the S-M. Finally we present some 
conclusions in Sec. IIXI 



II. SELF-MODULATION OF LEPTON AND 
PROTON BUNCHES 

The growth of the S-M instability for a long (k p cr z ^> 
kpz) bunch at position £ = z — ct along the bunch with 
peak density n^, after a propagation distance z in the 
plasma, is given by e G , where [T5] : 
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The expression for G from Ref. [15] is strictly valid for a 
flat top cylindrically symmetric narrow bunch (k p cr r <C 1) 
with matched emittance. For our parameters (k p cr r ~ 
0.9) there are almost no differences between the growth 
rates from Ref. [15] and Ref. [16], which is valid also for 
broad bunches with k p <j r > 1. 

Considering £ = and the same propagation dis- 
tance z, Eq. evaluated with the parameters of Ta- 
ble [I] shows that the growth rate in the lepton case is 
approximately 20 times larger than in the proton bunch 
case. The instability is therefore expected to develop 
over much shorter distances in the lepton case. Numeri- 
cal simulations support this evaluation, showing cm-scale 
evolution distances with leptons (see below) rather than 
meter-scale distances with protons [17]. A S-M experi- 
ment could therefore be performed using the long lepton 
bunches of early single bunch PWFA experiments [18] 
(ex 500 jam) and high-density plasmas currently avail- 
able for two ultra-short bunches PWFA experiments [9 
(oc 10 1T cm -3 ). With these parameters many of the S- 
M physics could be tested very soon over centimeter to 
meter-scale plasma lengths and with most of the diag- 
nostics also available (see below). In addition, seeding 
of the S-M instability could also be tested. Collimation 
techniques [19] that will be used to shape the short elec- 
tron bunch into a drive/witness train [9] can be used to 
create a sharp rising edge, shorter than the plasma wave- 
length, in the long bunch current distribution. A bunch 
with a sharp edge excites a low amplitude wake in a con- 
trolled manner making S-M more predictable. This is 
crucial to deterministically inject a witness bunch into 
the accelerating phase of the wakefields. 

Since short (« 20 jam) electron and positron bunches 
will also be available, acceleration of a witness bunch 
could be tested. Note that generation of wakefields and 
thus seeding of the instability could be accomplished by 
a short and intense laser pulse. Seeding by injecting the 
short electron or positron bunch in front of the long (op- 
posite charge) bunch could be explored. Moreover, with 
the pre-ionized plasma source that will be available [20] 
the sensitivity of the instability development and satu- 
ration to plasma density gradients could be determined 
[2T] , Plasma density variations larger than the inverse 
number of beamlets created along the bunch by the S- 
M instability are expected to detune the instability and 
decrease the wakefields amplitude. In addition, the de- 
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velopment of the instability along the plasma could be 
studied by varying the plasma length. 



Parameter 


PDPWFA 


PWFA 


<J r [/^m] 


200 


10 


<J r [c/Up] 


0.4 


0.9 


(j£ [cm] 


10 


5 x 1(T 2 




188 


45 


7o 


480 


4 x 10 4 


7o mc 2 /e[GeV] 450 


20.5 


-/Vpart 


11 x 10 iU 


2 x 10 1U 


n [cm -3 ] 


10 14 - 10 15 


2.3 x 10 17 


rib/ 'n 


2.0 x 10" 2 


io- 1 


L plasma [m] 


5 


1 


L plasma [c/cj p 


] 9 x 10 3 


9 x 10 4 


e n [mm • mrad] 3.83 


50 



TABLE I: Typical parameters for the long proton bunch 
experiments (PDPWFA) [15] and for the long electron and 
positron bunches experiments considered here (PWFA). The 
particle bunch length and transverse size are given by , and 
cr r , respectively. In addition, the plasma length is denoted by 
^plasma anc j e ^ ^ g ^ e normalized emittance. 



III. SIMULATION PARAMETERS 

We have performed 2D cylindrically symmetric and 3D 
numerical simulations using the fully relativistic, mas- 
sively parallel particle-in-cell code OSIRIS [22]. The 
2D (3D) simulations use a moving window that prop- 
agates at c, with resolutions k p Az = k p Ar = 0.0375 
(kpAz = kpAx = k p Ay = 0.2) for the longitudinal and 
transverse cell sizes, with 2x2 plasma and beam particles 
per cell. We use quadratic particle shapes, and current 



smoothing and compensation. Here k p = u p /c is the 



wake wave number, 
profile is: 
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for < £ - £ 

£o = cr^ a/27t . The bunch and plasma physical parameters 
are those of the PWFA case of Table [l| unless otherwise 
specified. The incoming bunch has zero initial energy 
spread. We do not present simulations for the proton 
bunch case, and refer to previous [15] results. 



IV. COMPETITION BETWEEN HOSING-LIKE 
AND SELF-MODULATION INSTABILITIES 

The 3D simulations show that the early propagation 
of the full length SLAC bunch is dominated by the com- 
petition between the S-M and the hosing-like instability. 




FIG. 1: Results of 3D OSIRIS simulations illustrating the 
logarithm of the bunch density after propagation in 7.5 cm 
of plasma. Figure (a) shows the competition between the 
hosing instability and S-M instability leading to breakup in 
the full length bunch case. Figure (b) shows that the hosing 
instability can be effectively mitigated by using a hard-cut 
bunch (cut in the middle of the bunch in this case) that seeds 
the S-M instability, enabling it to grow, saturate and further 
drive large amplitude wakefields over long plasma lengths. 
The bunches propagate towards the right as indicated by the 
arrows. Note that in these simulations rib /no = 5 x 10 -2 since 



part 



1 x 10 10 . 



The seed for the hosing in the simulations is the slight ra- 
dial asymmetries in the transverse distribution resulting 
from the initial macro-particle distribution. Actual par- 
ticle bunches may have incoming correlations (in space 
and momentum) and tilts that could seed the instability 
at a higher level. It is likely that these will be the domi- 
nant seed for hosing-like instabilities in experiments. The 
hosing-like instability leads to bunch break up shortly 
after the beginning of the plasma (Fig. [IJi), i.e., over 
a length scale comparable to that of the saturation of 
the S-M instability (~ 5 cm, see Sectior fVl| ). However, 
these simulations show that hosing can be mitigated by 
using bunches with a short rise time in comparison to 
the plasma wavelength X p = 2irc/ujp (Fig. JTJ)). This oc- 
curs because the short rise time of the beams provides 
an initial wakefield that can effectively seed the S-M in- 
stability so that it grows substantially before hosing oc- 
curs. Seeding mechanisms are important to stabilize the 
accelerating structures of the self-modulated PWFA [5] 
providing a controllable noise source and relative phase 
for the growth of the instability. 



V. ELECTRON AND POSITRON BUNCH 
SELF-MODULATION 

Now that we have established that seeding of the S-M 
instability can mitigate the growth of the hosing insta- 
bility, we use 2D cylindrically symmetric OSIRIS sim- 
ulations with assymetric bunches with rapid rise times 
to investigate the evolution of the wakefields over meter- 
long plasmas. We note that in the case of a bunch with 
a sharp rise time the initial wakefield amplitude is much 
higher than that of the full bunch. However, the num- 
ber of particles available to drive the instability is lower 
and the effective bunch length shorter. Therefore, the 
evolution of the self-modulated wakefields results from 
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FIG. 2: Cylindrically symmetric OSIRIS simulation results 
showing the electron bunch density after z—10 cm (a) and z 
= 1 m (b) and the positron bunch density after z=10 cm (c) 
and z = 1 m (d). The bunches propagate towards the right 
as indicated by the arrows. In both cases there is a hard-cut 
in the middle of the bunch at £ = 1.27 mm. Note that the 
density is shown in logarithmic scale. 



a balance between the initial amplitude of the seeding 
wakefields, and the corresponding number of bunch par- 
ticles available to drive the wakefield. This balance may 
be investigated in detail in experiments with long elec- 
tron and positron bunches and is briefly addressed in 
Section lYl! Bl for the electron bunch case. 

We start with the case of a half-cut bunch where the 
density profile is given by Eq. ^ with < £ < f - Fig- 
ure [2] shows the electron bunch density after a propaga- 
tion distance of (a) z = 10 cm and (b) z = 1 m (length 
comparable to that of the available plasma). These show 
that the S-M leads to the formation of several beam- 
lets (oc 17 in this case) that can be clearly identified in 
Fig. [2] a. Electrons initially propagating in (de) focusing 
regions of the wakefield are (de) focused, thus enhancing 
the initial (de) focusing plasma fields providing the feed- 
back for the transverse two-stream instability [5]. How- 
ever, just as for laser pulse self-modulation only half of 
the beamlet resides completely in a focusing region in the 
linear phase [7]. 

The positron bunch densities at z = 10 cm, and at z = 
1 m, are shown in Figs. [2] c and d. Figure [2] shows that 
the S-M instability associated with positron and electron 
bunches can differ due to nonlinear effects. The number 
of bunches is similar to that in the electron bunch case, 
however, there are already much less particles in the back 
of the bunch at z = 10 cm. Note that at this distance in 
the plasma the non-linear regime is already reached (see 
next Section). After the propagation of z = 1 m, the back 
of the positron bunch is strongly defocused, and the total 
number of positrons available to drive the wakefields is 
much lower than in the electron bunch case. Simulations 
show that after 1 m in the plasma 20% of the initial 
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FIG. 3: Maximum accelerating field E acce \ at r = a r (z = 
0)/2 (2D simulations) driven by the electron (blue line) and 
positron (red line) bunch along the plasma. The inset shows 
the logarithm of the fields along the first few centimeters, 
where the electric field grows exponentially in the linear 
regime of the instability and then saturates. The black line in 
the inset is the theoretical result for the growth rate (Eq. (fil)). 



number of positrons and 50% of the number of electrons 
are still driving the wake. Defocused particles typically 
leave the plasma wave at a relatively small angle, lower 
than 5 mrad. 



VI. WAKEFIELD EXCITATION 

A. Longitudinal wakefields 

Figure [3] illustrates the amplitude of the peak accel- 
erating fields as a function of the propagation distance 
in the plasma. It shows that the initial stage of the in- 
stability is similar for electron and positron bunches, in 
agreement with theoretical results [5j [15j [16] since the 
instability is in its linear growth phase. However, they 
differ once the wakefield amplitude reaches « 20 GV/m 
(^~5 cm), corresponding to the excitation of non-linear 
wakefields close to the blowout regime. For z > 5 cm the 
accelerating gradient driven by the positron bunch de- 
creases to nearly one-half of that of the electron bunch. 
Comparison between Eq. ([I]) (G at £ = cr^)and the simu- 
lation results is shown in the inset of Fig. (|3|, and we find 
very good agreement between theory and simulations for 
z < 4 cm, i.e. in the linear regime. 



B. Transverse wakefields 

The reason for the evolution of the accelerating field 
can be further explored by investigating the transverse 
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FIG. 4: Radial focusing (E r —cBe) field (black line, evaluated 
at r = a r (z = 0)/2, 2D simulation) and bunch charge density 
(red line) at z — 2.5 cm for (a) the electron and (b) the 
positron bunch. Only the £ = 0.55 mm to £ = 0.85 mm 
fractions of the bunches are plotted. The bunch propagates 
towards the right, as indicated by the arrows. 



radial electric field. Figure [4] shows the lineout of the 
radial plasma focusing force E r — cBq (evaluated at 
r = cr r o/2) where a r o is the initial spot-size, superim- 
posed with the lineout of the self-modulated electron 
(Fig. |4|l) and positron (Fig. ^p) bunch densities. These 
fields do not vary sinusoidally along z, as expected in the 
linear regime, even after such a short propagation dis- 
tance, confirming that non-linear, large amplitude wake- 
fields are excited. The nonlinear regime of the plasma 
wakefield is quickly reached in the seeded case consid- 
ered here because the initial wakefield amplitude (the 
seed) already reaches a few percents of the cold plasma 
wavebreaking field amplitude. In both cases, the leading 
edge of each nonlinearly self-modulated beamlet propa- 
gates in regions of defocusing fields (E r — cBq < for 
electrons, E r — cBq > for positrons), leading to the 
gradual erosion of the head of each beamlet (cf. Fig. [2)3). 
As a result the accelerating wakefield gradually decreases 
as it propagates (increasing z), which is consistent with 
Fig. [3] In the positron bunch scenario (Fig. |4Jd), as non- 
linear wakes are excited, the fields that are focusing for 
positrons (E r — cBq < for r > 0) are limited to the 
regions where the plasma electrons cross the axis in the 
back of each plasma period. As a result the beamlets 
are shorter in the positron case and more positrons are 
defocused. Hence the accelerating field is lower than in 
the electron bunch case. 



C. Energy gain and loss 



m 



The electron bunch energy spectrum at z = 1 
(Fig. |5| shows that the initial energy varies by up to 
15 GeV (loss at the 1%/GeV level of the peak value). 
For positrons, 8 GeV energy loss at the 1%/GeV level 
occurs. In both cases energy gain on the order of 5 GeV 
is observed. For the electron bunch case some beam par- 
ticles gain up to ~ 20 GeV. The broad continuous energy 
spectra are the result of the bunch particles sampling all 
the phases of the wakefield. 

It is interesting to note that as the blowout regime 
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FIG. 5: Electron (blue line) and positron (red line) bunch 
energy spectra after propagation of z — 1 m in the plasma. 
Both spectra are normalized to their peak value. 



is reached, the peak wakefield amplitudes driven by 
the long, self-modulated bunches are comparable to 
that obtained in recent PWFA experiments using com- 
pressed bunches of length approximately equal to that 
of the beamlets, but containing the full charge (« 
50 GV/m) [23U25]. This shows the effectiveness of the 
self-modulated bunch in driving large amplitude accel- 
erating fields through the resonant excitation of plasma 
wakefields. Uncompressed lepton bunches are more eas- 
ily produced in radio-frequency linacs and using a seeded 
self-modulated PWFA (rapid rise) may reduce the com- 
plexity of the drive bunch linac. However, since only 
about half the particles participate in the wake excitation 
and the wakefield amplitude varies along the bunch the 
energy transfer efficiency will be proportionally smaller. 



VII. EVOLUTION OF THE WAKEFIELDS 
ALONG THE PLASMA 

A. Phase velocity 

With available parameters (See Table[j]), the PDPWFA 
would operate in the linear regime [26] of the PWFA. Un- 
like the PDPWFA, the self-modulated PWFA can lead to 
the generation of non-linear wakes close to the blowout 
regime [13] . Because in the blowout regime the non-linear 
plasma wavelength differs from that of the linear regime, 
the phase velocity of the self-modulated wakefields suffers 
from additional spatial and temporal variations. How- 
ever, these issues can be overcome by splitting the S-M 
and the acceleration plasma sections. In an accelerator 
scheme, the self- modulation of the drive bunch may oc- 
cur in a short S-M plasma section, where the wake phase 
velocity is slow whe n com pared to that of the drive bunch 
[15J[T6] (see Section VII). The acceleration would then oc- 
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FIG. 6: Evolution of the accelerating wake structure (E acce \ 
field) within the bunch between £ = 0.55 mm and £ = 0.8 mm 
as a function of the propagation distance for the electron 
bunch case, (a) and for the positron bunch case (b). The 
dashed lines show the approximate trajectories of the peak 
accelerating field. 
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FIG. 7: (a) Maximum accelerating fields for the 3/4- 
cut bunch using electron (blue line) and positron (red line) 
bunches, (b) corresponding electron (blue line), and positron 
(red line) bunch spectra after 1 meter propagation. 



B. 1/2 versus 3/4 cut-bunch seeding 



cur in a following plasma section, where the wake phase 
velocity is close to that of the drive bunch. This scheme 
would allow for the external injection and trapping of a 
low energy witness bunch followed by its acceleration in 
the relativistic wakefields. Note that simulations show 
that external injection at an angle with respect to the 
drive bunch trajectory may lead to a narrow final energy 
spread even when the injected bunch is long compared to 
the plasma wavelength [15] . 

The qualitative evolution of the wakefield phase veloc- 
ity of lepton bunches can be examined through Fig. [6] 
It illustrates the evolution of the axial accelerating field 
amplitude driven by the electron (Figj6] a) and by the 
positron (Figj6] b) bunch. The phase velocity of the 
plasma waves is deduced from the slope of the lon- 
gitudinal electric field (E z ) in the z — £-plane. The slope 
is given by || = -^r c - Thus, a more (less) vertical slope 
corresponds to a faster (slower) phase velocity. Two re- 
gions can be identified on Fig. [6j First, before the sat- 
uration of the S-M instability, for z < 5 cm the wake 
phase velocity is significantly lower than the velocity of 
the beam ~ c. This is the region corresponding to 
the S-M plasma section. Second, after the saturation of 
the S-M instability, for z > 5 cm, the wake phase veloc- 
ity is close to ~ c and dephasing is not an issue [T5] . 
This is the acceleration plasma section. Figure [6] shows 
that in both regions is larger with the electron than 
with the positron bunch, further indicating that electron 
bunches are advantageous for the acceleration of an exter- 
nally injected witness particle bunch. Simulations show 
that these differences in self-modulated wakefields phase 
velocity are due to the smaller focusing field regions for 
positrons than for electrons in the blowout regime. This 
leads to faster erosion of the positron beamlets head 
(slower phase velocity) and lower overall wakefield am- 
plitudes. 



We also performed simulations to show how the results 
change when the location of the rapid rise was earlier in 
the bunch. We considered a bunch with a hard cut den- 
sity profile given by Eq. ^ with < f - £ < 3/2a^ V2n 
(3/4-cut bunch), rather than with < £ — £o < cr^ /V2tt 
(half-cut bunch), since this can be tested at SLAC 
FACET. The simulation results shown in Fig. [7^i reveal 
that the 3/4-cut bunch provides larger accelerating gra- 
dients in comparison to the half-cut bunch because it car- 
ries more charge and is longer, despite the lower initial 
amplitude of the seeding compared to the half-cut case. 
As a result, the corresponding particle energy spectra of 
Fig. [7]3 show energy variations in excess of 20 GeV at 
the 1%/GeV level for electrons, and 6 GeV for positrons 
(compare with Figs. [3] and [5|. Note however that the 
longer 3/4-cut bunch is more sensitive to the hosing-like 
instability than the 1/2-cut bunch. Its propagation sta- 
bility will be studied using 3D simulations in a future 
publication. 



VIII. OBSERVATION OF THE S-M 
INSTABILITY 

Since simulations indicate that the S-M and the ex- 
citation of wakefields reach the non-linear regime with 
available electron and positron bunches, we briefly ex- 
plore how S-M and its effects on the bunches could be 
diagnosed in an experiment. The S-M could be mea- 
sured with coherent radiation detection techniques. With 
the parameters considered here the period of the S-M 
is in the 70 jam range. Therefore, the coherent transi- 
tion radiation (CTR) emitted by the self- modulated train 
when traversing a thin metallic foil is in the 4 THz fre- 
quency range and below. Emission of CTR at frequen- 
cies > 1 THz when the bunch travels in the plasma as 
opposed to < 1 THz without plasma (a z = 500 /xm) 
would be a clear indication that S-M occurred. Such an 
integrated CTR energy measurement has already been 
successfully used to monitor the relative length of ultra- 



7 



short electron bunches on a shot-to-shot basis [25] . CTR 
auto-correlation measurements could also be used to de- 
termine the S-M period in an average sense, assuming 
that the development of the instability is sufficiently re- 
producible. Information about the bunch modulation 
depth and beamlets actual length could in principle be 
obtained from the spectrum of the coherent radiation. 
Fourier transforms of the electron and positron bunch 
longitudinal density profiles after 1 m of plasma are 
shown in Fig. [8j These could be obtained from single 
shot Smith-Purcell diagnostic [27, and show clear peaks 
at the plasma wavelength (k = k p in Fig. M). Peaks at 
the second harmonic (k = 2k p ) are also visible in Fig. |8j 
especially in the positron bunch case, characterized by 
shorter beamlets (see Fig. [2] and Fig. [4]). Unless the in- 
stability is seeded, single shot diagnostic are necessary 
to capture the instability-aspect of the self-modulation. 
Note however that in the case of ultra-relativistic bunches 
considered here collimation of the defocused particles to 
yield a real bunch current modulation at the diagnos- 
tic will be necessary. This is due to the fact that with 
high-energy particle the coherent radiation is essentially 
emitted in the forward direction, significantly decreasing 
the dependency of the emitted energy or of the emitted 
spectrum on the transverse size modulation. 

In addition, energy changes of the bunch particles 
could be easily measured since the values for energy 
gain and loss obtained from simulations are well within 
the energy measurement limits of previous single, short 
bunch PWFA experiments [23-25 j. Finally the focus- 
ing/defocusing effect of the S-M instability on the bunch 
particles could also be detected by observing the bunch 
transverse size, a short distance downstream from the 
plasma exit, using for example optical transition radia- 
tion imaging [23j [28] . 



by using full beams, and beams with rapid rise times that 
are predicted to mitigate its appearance. This is also con- 
firmed by numerical simulations. Other mechanisms that 
can lead to the suppression of the S-M instability, includ- 
ing the motion of the background plasma ions [29] , could 
also be investigated. The possible experiments present 
outstanding opportunities to unravel the physical mech- 
anisms associated with the propagation of long particle 
beams in dense plasmas. 




k[k P ] 

FIG. 8: Absolute value of the fast Fourier transform of the 
electron (blue line) and positron (red line) bunch on-axis den- 
sity after 1 m of plasma to be compared with that of the un- 
modulated incoming bunch (black line). The peak at k = k p 
corresponds to the S-M of the bunches at the plasma wave- 
length. The peak at the first harmonic (k — 2k p ) is visible in 
the positron bunch case, for which the beamlets are shorter. 
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